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Abstract —Nonlinear analysis of the phase-locked loop 
(PLL) based circuits is a challenging task, thus in 
modern engineering literature simplified mathematical 
models and simulation are widely used for their study. 
In this work the limitations of numerical approach is 
discussed and it is shown that, e.g. hidden oscillations 
may not be found by simulation. Corresponding exam¬ 
ples in SPICE and MatLab, which may lead to wrong 
conclusions concerning the operability of PLL-based 
circuits, are presented. 


I. Introduction 

The phase-locked loop based circuits (PLL) are widely 
used nowadays in various applications. PLL is essentially 
a nonlinear control system and its rigorous analytical 
analysis is a challenging task. Thus, in practice, simulation 
is widely used for the study of PLL-based circuits (see, e.g. 
0 "@)- At the same time, simulation of nonlinear control 
system may lead to wrong conclusions, e.g. recent work 5] 
notes that stability in simulations may not imply stability 
of the physical control system, thus stronger theoretical 
understanding is required. 

In this work the two-phase PLL is studied and corre¬ 
sponding examples, where simulation leads to unreliable 
results, is demonstrated in SPICE and MatLab. 

II. PLL OPERATION 

Typical analog PLL consists of the following elements: 
a voltage-controlled oscillator (VCO), a linear low-pass 
filter (LPF), a reference oscillator (REF), and an analog 
multiplier ® used as the phase detector (PD). The phase 
detector compares the phase of VCO signal against the 
phase of reference signal; the output of the PD (error 
voltage) is proportional to the phase difference between its 
two inputs. Then the error voltage is filtered by the loop 
filter (LPF). The output of the filter is fed to the control 
input of the VCO, which adjusts the frequency and phase 
to synchronize with the reference signal. 

Consider a signal space model of the classical analog 
PLL with a multiplier as a phase detector (see Fig. [lj. 
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Fig. 1: Operation of classical phase-locked loop for sinusoidal 
signals 


Suppose that both waveforms of VCO and the reference 
oscillator signals are sinusoidal [|] (see Fig. [T]) . The low- 
pass filter passes low-frequency signal sin(0i(ij— 6L(t)) and 
attenuates high-frequency signal sin(#i(t) + doft)). 

The averaging under certain conditions |6|, 7 , [9|- 
11 and approximation tp(t) ss sin(0i(i) — 62 ( 1 )) allow one 
to proceed from the analysis of the signal space model 
to the study of PLL model in the signal’s phase space. 
Rigorous consideration of this point is often omitted (see, 
e.g. classical books 12 p.l2,pl5-17] •.ESI, p/7]) while it may 
lead to unreliable results (see, e.g. [l4], [l5]). 

One of the approaches to avoid this problem is the use of 
two-phase modifications of PLL, which does not have high- 
frequency oscillations at the output of the phase detector 
16 . 


III. Two-phase PLL 
Consider two-phase PLL model in Fig. [2] 



Fig. 2: Two-phase PLL 


Here a carrier is cos(0i(i)) with 9i(t) as a phase and 
the output of Hilbert block is sin(0i(f)). The VCO gen¬ 
erates oscillations sin($ 2 (f)) and cos(6L(t)) with & 2 {t) as 
a phase. Fig. [3] shows the structure of phase detector 
(complex multiplier). The phase detector consists of two 

1 Other waveforms can be similarly considered 
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si n (0-, (t)) 



analog multipliers and analog subtractor. The output of 

PD is ip(t) = sin(0i(f))cos(02(i)) — cos(0i(f))sin(02(i)) = 
sin(0i(t) — 02 (i)) In this case there is no high-frequency 
component at the output of phase detector. It is reasonable 
to introduce a phase detector gain \ (e.g. to consider 
additional loop filter gain equal to 0.5) to make it the same 
as a classic PLL phase detector characteristic: 

<p(t) = ^sin(6>i(t)-0 2 (t)). (1) 

Consider a loop filter with the transfer function H(s). 
The relation between input ip(t) and output g(t ) of the 
loop filter is as follows 

x = Ax + bip(t), g(t) = c* x + h<p(t ), 

H(s) = c*{A-siy l b~h. 

The control signal g(t) is used to adjust the VCO phase 
to the phase of the input carrier signal: 

^2 (r)dT = Wf ree t + L / g{r)dr, (3) 

J o 

where o4 ree is the VCO free-running frequency (i.e. for 
g(t) = 0) and L the VCO gain. 

Next examples show the importance of analytical meth¬ 
ods for investigation of PLL stability. It is shown that the 
use of default simulation parameters for the study of two- 
phase PLL in MatLab and SIMULINK can lead to wrong 
conclusions concerning the operability of the loop, e.g. 
the pull-in (or capture) range (see discussion of rigorous 
definitions in [17 , [l8]). 



been an attractive theoretical problem. This problem has not 
yet been provided with a satisfactory solution". Below we 
demonstrate that in this case a numerical simulation may 
give wrong estimates and should be used very carefully. 

In Fig.[4]we use the block Loop filter to take into account 
the initial filter state a:(0); the initial phase error 0^(0) 
can be taken into account by the property initial data 
of the Intergator blocks]^] Note that the corresponding 
initial states in SPICE (e.g. capacitor’s initial charge and 
inductor’s initial currents) are zero by default but can be 
changed manually. 

In Fig. [5] the two-phase PLL model simulated with 
relative tolerance set to “le-3” or smaller does not ac¬ 
quire lock (black color), but the PLL model in signal’s 
phase space simulated in MatLab Simulink with standard 
parameters (a relative tolerance set to “auto”) acquires 
lock (red color). Here the input signal frequency is 10000, 
the VCO free-running frequency u4 ree = 10000 —178.9, the 
VCO input gain is L = 500, the initial state of loop filter is 
xq = 0.1318rl and the initial phase difference is 9/± (0) = 0. 



Fig. 5: Simulation of two-phase PLL. Filter output g(t) for the 
initial data xq = 0.1318,$a( 0) =0 obtained for default “auto” 
relative tolerance (red) — acquires lock, relative tolerance set 
to “le-3”(green) — does not acquire lock. 


V. Simulation in SPICE 


IV. Simulation in MatLab 

Consider a passive lead-lag loop filter with the transfer 
function H(s) = 1+ ^^ T2 ) , n = 0.0448, r 2 = 0.0185 and 


the corresponding parameters A = — 




c = 


h = T2 

T\ +T2 ’ n +T2 ‘ 


T1+T2 ’ T1+T2’ 

The model of two-phase PLL in 


' X I ' Z L I Z I I 

MatLab is shown in Fig. 141 (see more detailed description 
of simulating PLL based circuits in MatLab Simulink in 

H-P). 

For the 


case of the passive lead-lag filter a recent work 
22 p. 123] notes that “the determination of the width of 
the capture range together with the interpretation of the 
capture effect in the second order type-I loops have always 


In this section the previous example is reconstructed 
in SIMetrix, which is one of the commercial versions of 
SPICE. 

Consider SIMetrix model of two-phase PLL shown in 
Fig- ID, The input signal and the output of Hilbert block 
in Fig. [2] are modeled by sinusoidal voltage sources VI (a 

2 Following the classical consideration [l2| p.17, eq.2.20] |l3| p.41, 

eq.4-26], where the filter’s initial data is omitted, the filter is often 
represented in MatLab Simulink as the block Transfer Fen with zero 
initial state (see, e.g. [23] — [2T| ). It is also related to the fact that the 
transfer function (from ip to g) of system is defined by the Laplace 
transformation for zero initial data ir(0) = 0. 

3 Almost each initial state from the interval [1,2] gives similar 
results. 
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Trigonometric vco free-running frequencyl 


Function5 

Fig. 4: Model of two-phase PLL in MatLab Simulink 



frequency parameter is 1.5915494fc) and V2 (a frequency 
parameter is 1.5915494/c and a phase is 90) (sin_input and 
cos_input). A complex multiplier in Fig. Ons modeled as 
two arbitrary sources ARBI and ARB2 with definitions set 
to V(N1)*V(N2). To subtract the output signals of mul¬ 
tipliers, Voltage Controlled Voltage Source (E3) is used. 
Phase detector gain (E5) is equal to Loop Filter in Fig. [2^ 
is modeled as a passive lead-lag filter with resistor R27 
capacitor C2, and resistor Rl. The input gain of VCO (E6) 
is equal to —500. VCO self frequencj|2](DC Voltage Source 
V3) is set to 9.8211fc. Voltage Controlled Voltage Source 
E2 summarizes a VCO self frequency and a control signal 
from E6. Resistor Rlbl (10u), capacitor Cl (5G), and am¬ 
plifier El(50fc) form an integrator)] The VCO waveforms 
are defined by arbitrary blocks ARB3 (with the function 


4 Zero input response (ZIR) frequency 

5 The same results could qualitatively be obtained using 10 K for 
the resistor and 5 Farad for the capacitor (keeping RC constant) 


sin(V r (A r l))) and ARB4 (with the function cos(V(iVl)))). 
Netlist for the model, generated by SIMetrix, is as follows: 


1 

2 

3 

4 

5 

6 
7 


8 

9 

10 

11 


12 


*#SIMETRIX 

VI sin_Input 0 0 Sine(0 1 1.5915494k 0 
c i 0 ) 

V2 cos_input 0 0 Sine(0 1 1.5915494k 
^ -157.03518u 0) 

V3 vco_frequency 0 9.8211k 
Rl C2_N 0 1.85k 
R2 filter_out PD_output 4.48k 
X$ARB1 sin_Input vco_cos_output 
'->• ARB1_0UT $$arbsourceARBl 
c ->- pinnames: N1 N2 OUT 
.subckt $$arbsourceARBl N1 N2 OUT 
B1 OUT 0 V=V(N1)*V(N2) 

. ends 

X$ARB2 cos_input vco_sin_output E3_CN 
$$arbsour ceARB2 pinnames: N1 N2 
OUT 

.subckt $$arbsourceARB2 N1 N2 OUT 
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13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 


32 

33 

34 


B1 OUT 0 V=V(N1)*V(N2) 

. ends 

X$ARB3 integrator_out vco_sin_output 

°4 $$arbsourceARB3 pinnames: N1 OUT 
.subckt $$arbsourceARB3 N1 OUT 
B1 OUT 0 V=sin(V(Nl)) 

. ends 

X$ARB4 integrator_out vco_cos_output 

°4 $$arbsourceARB4 pinnames: N1 OUT 
.subckt $$arbsourceARB4 N1 OUT 
B1 OUT 0 V=cos(V(Nl)) 

. ends 

El integrator_out 0 E1_CP 0 50k 
E2 integrator_in 0 vco_frequency E2_CN 

^4 1 

Cl E1_CP 0 5G 

C2 filter_out C2_N lOu IC=185m 

^4 BRANCH={IF(ANALYSIS=2,1 ,0)} 

E3 E3_P 0 ARB1_0UT E3_CN 1 
E5 PD_output 0 E3_P 0 500m 
E6 E2_CN 0 filter_out 0 -500 
Rlbl integrator_in E1_CP lOu 
.GRAPH filter_out curveLabel= 

°4 filter_out nowarn=true ylog=auto 
°4 xlog=auto disabled=false 
.TRAN 0 5 0 lm UIC 
.OPTIONS minTimeStep=lm 
+ tnom=27 


In Fig. [7] are shown simulation results in SPICE, which 
are close to the simulation results in MatLab Simulink (see 
Fig. [§ . For default simulation parameters in SIMetrix two- 
phase PLL synchronizes to the reference signal (red line). 
However, if we choose smaller simulation step (lm), the 
simulation reveals an oscillation (green line). 


“i 



Time/Secs 


ISecs/div 


Fig. 7: Hidden oscillations in SPICE 


VI. Mathematical reasoning 

Two-phase PLL is described by ([!]), p|, and Q, which 
form the following system of differentialequations 

x = Ax+ - sin(0/^), 

9 A =w A -Lc*x-^sm(6 A ), ^ 

9 A (t) = 9 1 (t)-9 2 (t), uj a = uji ojf ree 

For a lead-lag filter, described by the transfer function 
H(s) = l+gfti+r^) ) system Q takes the form 


-1 


+ ( 1 - 


T2 ) a sin(0 A ), 


ri+r 2 " ' v ri + r 2 2 

9 a =uj a -L - x -—— — sin(0 A ). 

ri + r 2 n + r 2 2 


(5) 


The equilibrium points of ^ are defined by the following 
relations: 


%eq = ySin(6» A ), sin(9 eq ) = 2-^-. 


( 6 ) 


For n = 0.0448, L = 500, and uj a = 178.9 we get 


x e „ = 0.016, 

9 (7) 

9 eq = (—l) fe 0.7975 + 7rfc, k G N. w 

Consider now a phase portrait (where the system’s 
evolving state over time traces a trajectory (x(tY9 A (t))), 
corresponding to signal’s phase model (see Fig. [8]). 



Fig. 8: Phase portrait of the classical PLL with stable and 
unstable periodic trajectories 


The solid blue line in Fig. [^corresponds to the trajectory 
with the loop filter initial state cc(0) = 0.005 and the VCO 
phase shift 0 rad. This line tends to the periodic trajectory, 
therefore it will not acquire lock. All the trajectories under 
the blue line (see, e.g., a green trajectory with the initial 
state x(0) = 0) also tend to the same periodic trajectory. 

The solid red line corresponds to the trajectory with 
the loop filter initial state 0.00555 and the VCO initial 
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— real trajectories simulation 



Fig. 9: Phase portrait of the classical PLL with stable and 
unstable periodic trajectories 


another stable solution may depend on the initial data and 
integration step. 


VII. Conclusion 


The considered example is a motivation for the use of 
rigorous analytical methods for the analysis of nonlinear 
PLL models. Rigorous study of the above effect can be 
done by Andronov’s point transformation method and 
phase plane analysis. Correspon ding bifurcation diagram 


were given in 49 (see, also 11 , 28 ). 


For Costas loop models, similar effect is shown in 19 


20 , 50 


phase 0. This trajectory lies above the unstable periodic 
trajectory and tends to a stable equilibrium. In this case 
PLL acquires lock. 

All the trajectories between stable and unstable periodic 
trajectories tend to the stable one (see, e.g., a solid green 
line). Therefore, if the gap between stable and unstable 
periodic trajectories is smaller than the discretization 
step, the numerical procedure may slip through the stable 
trajectory. The case corresponds to the close coexisting 
attractors and_the bifurcation of birth of semistable tra¬ 
jectory 28 , 29 . In this case numerical methods are 


limited by the errors on account of the linear multistep 
integration methods (see [30 , 31 ). As noted in [32], low- 
order methods introduce a relatively large warping error 
that, in some cases, could lead to corrupted solutions (i.e., 
solutions that are wrong even from a qualitative point 
of view). This example demonstrate also the difficulties 
of numerical search of so-called hidden oscillations, whose 
basin of attraction does not overlap with the neighborhood 
of an equilibrium point, and thus may be difficult to 
find numerical!}]^] In this case the observation of one or 


6 An oscillation in a dynamical system can be easily localized 
numerically if the initial data from its open neighborhood lead to 
long-time behavior that approaches the oscillation. From a compu¬ 
tational point of view, on account of the simplicity of finding the 
basin of attraction in the phase spa ce, i t is natur al to suggest the 
following classification of attractors [28] , [33] — [36] : An attractor is 
called a hidden attractor if its basin of attraction does not intersect 
small neighborhoods of equilibria, otherwise it is called a self-excited 
attractor. 

For a self-excited attractor its basin of attraction is connected 
with an unstable equilibrium. Therefore, self-excited attractors can 
be localized numerically by the standard computational procedure in 
which after a transient process a trajectory, started from a point 
of unstable manifold in a neighborhood of unstable equilibrium, is 
attracted to the state of oscillation and traces it. Thus self-excited 
attractors can be easily visualized. 

In contrast, for a hidden attractor its basin of attraction is not con¬ 
nected with unstable equilibria. For example, hidden attractors can 
be attractors in the systems with no equilibria or with only one stable 
equilibrium (a special case of multistable systems and coexistence of 
attractors). Recent examples of hidden attractors can be found in The 
European Physical Journal Sp ecia l Top ics: Multistability: Uncovering 
Hidden Attractors , 2015 (see [37] — [48] ). 
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